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Abstract 

Refractory  diboride  with  silicon  carbide  additive  has  a  unique  oxide  scale  structure  with  two  condensed 
oxide  phases  (solid  +  liquid),  and  demonstrates  oxidation  resistance  superior  to  either  monolithic  diboride 
or  silicon  carbide.  We  rationalize  that  this  is  because  the  silica-rich  liquid  phase  can  retreat  outward  to 
remove  the  high  SiO  gas  volatility  region,  while  still  holding  onto  the  zirconia  skeleton  mechanically  by 
capillary  forces,  to  form  a  ’’solid  pillars,  liquid  roof’  scale  architecture  and  maintain  barrier  function. 
Basic  assessment  of  the  oxygen  carriers  in  the  borosilicate  liquid  in  oxygen-rich  condition  is  performed 
based  on  first-principles  calculations.  It  is  estimated  from  entropy  and  mobility  arguments  that  above 
a  critical  temperature  Tc  ~  1500°C,  the  dominant  oxygen  carriers  should  be  network  defects,  such  as 
peroxyl  linkage  or  oxygen  deficient  centers,  instead  of  molecular  as  in  the  Deal-Grove  model.  These 
network  defects  will  lead  to  sub-linear  dependence  of  the  oxidation  rate  with  external  oxygen  partial 
pressure. 

Analysis  of  Microstructural  Effect  in  Oxidation 

Oxidation  of  ZrB2+SiC  produces  a  complex  scale  structure  containing  at  least  two  condensed  phases: 
Zr02  (c)  which  has  highly  porous  skeleton  with  percolating  holes,  and  a  silica-rich  liquid  phase  that  wets 
the  skeleton.  Gas  species  of  the  greatest  interest  are  B2O3,  SiO  and  CO,  although  BO,  B2O2,  B2O  etc. 
are  also  present  and  can  play  important  roles  in  mass  transport.  SiO  could  evolve  by  for  instance:  2Si02 
(1)  — >  2SiO  (g)  +  O2,  which  is  a  key  reaction  that  has  been  used  in  constructing  volatility  diagrams. 
When  in  contact  with  Si02  (1),  with  decreasing  oxygen  chemical  potential  or  the  equivalent  partial 
pressure,  SiO  will  have  higher  equilibrium  vapor  pressure.  Volatility  diagram  of  the  ZrB2+SiC  system 
indicates  that  when  T  >  Tb  —  1775°C,  the  peak  equilibrium  vapor  pressure  of  SiO  inside  the  scale  could 
exceed  1  atm,  which  would  then  induce  a  boiling  transition  (gas  bubbles  can  nucleate  and  grow  inside 
the  liquid).  This  violently  disrupts  the  Si02  (1)  scale  in  the  case  of  oxidizing  monolithic  SiC.  However, 
the  scales  of  ZrB2+SiC  and  HfB2+SiC  appear  to  be  much  more  tolerant  of  such  a  boiling  transition. 
It  is  precisely  in  this  T  >  Tb  regime  that  ZrB2-fSiC  and  HfB2+SiC  demonstrate  oxidation  resistance 
superior  to  monolithic  SiC. 

We  suggest  a  ’’dynamic  view”  and  a  ’’steady-state  view”  of  why  ZrB2+SiC  is  superior  to  monolithic 
SiC.  The  two  views  are  inherently  consistent.  Imagine  a  ZrB2+SiC  specimen  is  gradually  being  heated 
up  in  an  oxygen-rich  environment  like  normal  air  (Po2=0.2  atm),  initially  from  T  <3C  Tb-  At  such  high 
ambient  Po2,  a  protective  Si02  (1)  film  will  condense  on  top  from  the  very  beginning,  that  wets  the  ZrC>2 
(c)  skeleton,  leaving  no  voids  at  the  base.  At  T  <  Tb,  monolithic  SiC  in  fact  resist  oxidation  better 
than  ZrB2+SiC.  Also  at  T  <  Tb,  the  Psio  branch  of  the  volatility  diagram  in  contact  with  SiC>2  (1)  has 
a  formal  thermodynamic  definition  but  is  not  physically  realizable  as  pure  SiO  gas  bubbles,  because  any 
SiO  gas  bubble  anywhere  will  be  crushed  by  the  dual  forces  of  surface  tension  and  hydrostatic  pressure, 
which  we  take  to  be  1  atm  inside  Si02  (1).  However,  as  the  temperature  is  brought  up  to  T  >  Tb,  a 
sharp  transition  happens  inside  Si02  (1).  Now  SiO  gas  bubbles  can  nucleate  at  the  base,  with  PsiO-latm 
>  0  working  against  the  surface  tension.  The  dynamic  view  examines  how  the  gas  bubbles  grow  and 
coalesce,  paying  attention  to  the  role  of  the  Zr02  (c)  skeleton. 
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The  ZrC>2  (c)  skeleton  will  regulate  gas  bubble  dynamics.  Unlike  unconstrained  growth  inside  a  com¬ 
pletely  liquid  scale,  gas  bubbles  (SiO,  CO,  B2O3,  BO,  O2  etc.  mixture)  in  a  semi-solid  porous  scale 
arc  forced  to  grow  into  long  fingers  (a  pressure  difference  of  the  order  atm  is  large  enough  to  displace 
a  liquid,  but  usually  not  enough  to  displace  a  solid).  Reaction  2Si02  (1)  — ►  2SiO  (g)  +  O2  could  then 
happen  on  one  end  of  the  gas  finger,  BO,  B2O3,  SiO,  CO  etc.  would  then  diffuse  along  the  gas  finger 
with  O2  diffusing  in  the  opposite  direction,  and  finally  when  Po2  gets  high  enough,  SiO  could  get  re¬ 
oxidized  to  form  Si02  (1)  on  the  other  end  of  the  finger,  and  B2O3  etc.  would  get  solvated  in  the  liquid 
and  continue  to  diffuse  up  the  scale.  This  is  equivalent  to  a  channeling  transfer  of  Si02  (1)  from  one 
end  of  the  finger  to  the  other,  which  is  mechanically  untenable  without  the  support  and  constraint  of 
the  Zr02  (c)  skeleton.  In  reality  the  pores  are  tortuous  instead  of  straight,  giving  the  effusing  SiO  (g) 
much  opportunity  to  react  with  O2  (g)  near  the  end,  and  the  Si02  (1)  product  collected  on  the  Zr02 
(c)  substrate.  The  skeleton  may  also  impart  significant  mechanical  integrity  to  the  scale  in  the  case  of 
bubble  outbreaks  (bubble  diameter  constrained  by  the  pore  diameter)  or  under  external  shear  flow,  since 
Si02  (1)  adheres  strongly  to  a  highly  porous  Zr02  (c)  skeleton  by  capillary  forces  across  a  large  contact 
area.  In  short,  the  dynamical  view  is  that  the  Zr02  (c)  skeleton  helps  to  collect  and  retain  the  silica-rich 
liquid,  playing  an  important  mechanical  role. 

There  is  also  a  ’’steady-state”  explanation,  beginning  with  the  interpretation  of  volatility  diagrams. 
Volatility  diagrams  represent  chemical  equilibria  when  assuming  some  volatile  gas  species  are  in  contact 
with  certain  condensed  phases  -  solid  or  liquid.  Solid  or  liquid  makes  a  difference  here,  because  a  liquid 
could  flow  in  or  flow  out,  easily  retreating  from  a  region  if  necessary.  If  a  certain  condensed  phase 
retreats,  then  an  originally  high  volatile  gas  pressure  -  assuming  the  condensed  phase  was  there  -  loses 
its  significance.  This  understanding  of  the  volatility  diagram  may  be  translated  into  the  following  rule: 
liquid  phases  will  retreat  from  the  region  where  some  volatile  gas  species,  were  they  in  contact,  will 
have  high  vapor  pressure  (approaching  hydrostatic  pressure  inside  the  liquid),  to  region  of  lower  vapor 
pressure  according  to  the  volatility  diagram;  after  the  retreat,  the  actual  vapor  pressure  of  the  volatile 
gas  species  will  be  automatically  lowered  than  what  the  volatility  diagram  has  originally  indicated  for 
the  evacuated  region. 

The  above  rule  comes  from  thermochemistry.  Applying  the  rule  to  oxidizing  monolithic  SiC  above  T# , 
we  see  that  the  only  thermochemically  sound  steady-state  arrangement  is  for  Si02  (1)  to  retreat  to  the 
low  PSiO  ~  high  Po2  region,  and  let  gas-phase  diffusion  take  over  in  the  intervening  gap,  where  there  will 
be  no  condensed  liquid  phase,  thus  shutting  down  the  high  gas  volatility.  Unfortunately,  although  this 
setup  is  thermochemically  and  diffusion-kinetically  sound,  it  is  obviously  mechanically  unstable.  The 
’’scale”  would  easily  shear  off  or  spall.  This  is  fundamentally  because  in  monolithic  SiC,  with  only  a 
single  condensed-phase  oxide  product,  which  is  a  liquid,  there  is  no  way  to  satisfy  both  thermochemical 
and  mechanical  stabilities  simultaneously  at  T  >  Tq. 

In  contrast,  when  oxidizing  ZrB2+SiC,  one  gets  two  condensed-phase  oxide  products.  Zr(>2  (c)  itself 
has  only  low  volatilities  of  ZrO  (g)  and  ZrC>2  (g)  Furthermore  it  is  a  solid.  So  it  does  not  retreat  from 
the  high  SiO  volatility  region,  maintaining  mechanical  connection  with  the  main  body.  The  silica-rich 
borosilicate  liquid  duly  retreats  from  the  high  SiO  volatility  region,  thereby  removing  the  high  SiO 
volatility  automatically.  This  staggered  placements  of  solid  and  liquid  phases  at  T  >  Tb,  with  internal 
gas  finger  diffusion  and  gas-solid  reactions,  is  both  thermochemically  and  mechanically  sound  if  the  Zr02 
(c)  skeleton  is  long  enough,  such  that  even  after  the  retreat  the  liquid  phase  can  still  hold  onto  the  solid 
by  capillary  forces,  to  have  a  ’’solid  pillars,  liquid  roof’  architecture  as  suggested  by  many  experiments. 
The  actual  oxidation  of  SiC  particles  at  T  >  Tb  then  occur  directly  by  SiC  (c)  +  O2  (g)  — 1 ►  SiO  (g) 
4-  CO  (g),  a  gas-solid  reaction  without  going  through  the  liquid  phase,  which  will  lead  to  a  porous 
“SiC-depleted”  substrate  layer  in  the  ZrB2+SiC.  The  Zr02  (c)  “solid  pillars”  will  also  grow  longer  at 
the  base  by  gas-solid  reaction,  with  O2  (g)  reactant  and  B*Oy  (g)  product.  On  the  other  end  of  the 
gas  finger,  we  will  have  the  reverse  reaction:  2SiO  (g)  +  O2  -»  2Si02  (1)  which  replenishes  the  ’’liquid 
roof’. 
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The  above  analysis  rationalizes  why  the  microstructure  of  ZrB2+SiC  is  important  for  its  oxidation 
resistance,  since  the  pore  sizes  of  ZrC>2  (c)  might  be  related  to  the  pre-oxidation  SiC  particle  sizes.  As 
Gasch  et  al.  mentioned,  “at  20  volume  percent  SiC,  if  the  SiC  particles  are  assumed  to  be  small  spheres 
randomly  distributed  throughout  the  HfB2  matrix,  the  amount  of  SiC  should  be  above  the  percolation 
threshold.  This  means  that  the  SiC  particles  form  a  network  that  is  interconnected  in  three  dimensions.” 
For  a  certain  fixed  volume  fraction,  smaller  pores  and  better  connectivity  inside  ZrC>2  (c)  could  enhance 
the  collection  and  retention  of  the  silica-rich  liquid.  Thus,  nanoscale  SiC  particles  might  improve  the 
oxidation  resistance  of  ZrB2+ SiC,  by  refining  the  microstructure  of  the  in  situ  formed  Zr02  (c)  skeleton1 . 

The  proposed  view  also  explains  why  70-80  vol%  of  the  composite  is  dedicated  to  ZrB2-  It  is  seen  that  a 
porous  Zr02  could  be  advantageous  for  the  oxidation  resistance  for  a  mechanical  reason  (not  a  diffusion 
kinetics  one),  if  there  is  also  a  liquid  oxide  product  phase  to  ’’collaborate”  with.  The  porous  skeleton 
needs  to  be  strong  enough  as  well  as  sufficiently  long,  to  have  enough  room  for  the  liquid  oxide  phase 
to  retreat  outward.  Otherwise,  the  two  phases  may  still  have  to  separate  (“high  volatility  blows  away 
the  liquid  roof”),  and  the  entire  system  would  lose  oxidation  protection.  One  cannot  have  too  much  SiC 
(and  thus  the  liquid)  and  not  enough  oxide  skeleton,  and  maintain  the  mechanical  and  thermochemical 
stabilities  of  the  “solid  pillars,  liquid  roof”  architecture.  Since  capillary  force  holds  the  solid  and  liquid 
phases  together,  microstructural  refinement  of  the  oxide  scale  will  lead  to  stronger  capillary  adhesion 
per  unit  volume,  which  could  lead  to  significant  improvement  of  the  overall  oxidation  resistance. 

Modeling  of  Oxygen  Transport  in  Silica-Rich  Liquid 

From  the  above  analysis,  we  see  that  if  Zr02  (c)  has  a  percolating-holes  microstructure,  the  borosilicate 
liquid  will  define  the  effective  barrier  against  oxygen,  irrespective  of  whether  fully  dense  ZrC>2  is  intrin¬ 
sically  better  barrier  (in  open-circuit  condition)  against  oxygen  or  not  than  the  borosilicate  liquid,  as 
the  liquid  occupies  both  serial  and  parallel  oxygen  transport  pathways  in  the  ’’solid  pillars,  liquid  roof’ 
architecture.  In  this  section  we  focus  on  the  atomic-level  events  that  govern  oxygen  transport  in  the 
borosilicate  liquid. 

Experimentally,  it  is  still  challenging  to  accurately  determine  the  composition  profile  of  the  borosilicate 
liquid  because  boron  is  a  light  element.  According  to  the  Hertz- Knudsen- Langmuir  equation,  the  net 
evaporation  flux  of  a  species  from  a  liquid  surface  is  proportional  the  difference  between  the  equilibrium 
vapor  pressure  and  the  actual  vapor  pressure  of  the  species  at  the  surface.  Since  pure  B2O3  (1)  has  much 
higher  equilibrium  vapor  pressure  than  pure  Si02  (1)  in  an  oxygen-rich  environment,  a  borosilicate  liquid 
facing  air  would  preferentially  evaporate  B2O3  instead  of  Si02-  Thus  the  borosilicate  liquid  should  be 
overall  silica  rich,  with  a  composition  gradient  that  is  B203-depleted  at  the  liquid-air  interface,  and 
B2C>3-€nriched  at  the  gas  finger-liquid  interface,  as  B2O3  (g)  and  other  boron-bearing  gas  species  are 
carried  along  with  SiO  (g)  in  the  gas  finger  and  get  absorbed  into  the  liquid.  Previously,  Bongiorno 
and  Pasquarello  have  studied  oxygen  transport  in  pure  silica  glass  using  a  multiscale  modeling  approach 
that  combines  high-level  quantum  mechanical  (DFT)  calculations  of  the  diffusing  oxygen  species  and 
local  energy  barriers,  with  kinetic  Monte  Carlo  sampling  of  connected  migration  pathways.  To  bound 
the  results,  we  decide  to  model  a  borosilicate  liquid  composition  of  equal  B2O3  and  SiC>2  proportions. 

One  major  challenge  in  modeling  any  glass  or  liquid  is  to  have  reliable  atomic  structures.  We  have 
adopted  the  structure  generation  approach  of  Van  Ginhoven,  Jonsson  and  Corrales,  which  was  shown 
to  reproduce  experimental  pair  distribution  functions  for  pure  silica.  The  approach  requires  a  classical 
interatomic  potential  to  perform  long-time  molecular  dynamics  (MD)  simulations  at  the  beginning, 
followed  by  further  DFT  optimizations.  To  generate  the  classical  potential,  we  adopt  the  van  Beest, 
Kramer  and  van  Santen  (BKS)  parameterization  for  Si-O  interactions,  but  fit  B-O  and  Si-B  interactions 
to  a  series  of  small  DFT  calculations  for  bulk  B2O3,  using  the  software  package  GULP.  Then,  starting 
from  random  positions  for  oxygen,  boron  and  silicon  atoms  in  the  supercell,  we  perform  a  sequence 
of  classical  MD  simulations  at  temperatures  6000,  5000,  4000,  3000,  2000  and  1000  K.  The  resulting 

1S.  S.  Hwang.  A.  L.  Vasiliev.  and  N.  P.  Padture,  Materials  Science  and  Engineering  A,  464  (2007)  216-224. 
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structures  were  then  used  as  input  geometries  for  further  DFT  calculations  (Vienna  ab-initio  simulation 
program  with  spin-polarized  PW91  functional,  projector  augmented  wave  method,  planewave  kinetic 
energy  cutoff  400  cV).  For  our  initial  studies,  a  cubic  supercell  containing  14  SiC>2  +  7  B2O3  formula 
units  is  used,  with  total  77  atoms.  The  average  density  is  2.3  g/cm3.  A  typical  liquid  structure  at  T  = 
2500°C,  after  further  equilibration  by  ab  initio  MD,  clearly  has  a  framework  structure  with  no  dangling 
bonds  (all  Si  are  4-fold  coordinated  to  O,  and  all  B  are  3- fold  coordinated  to  O). 

Based  on  these  atomic  structures,  we  have  studied  the  thermodynamic  stability  and  diffusion  kinetics 
of  solvated  oxygen  molecules  OJ  and  atomic  O*  in  borosilicate  liquid.  The  former  stay  inside  the  open 
cages  of  the  framework  and  do  not  interact  chemically  with  the  framework.  The  latter  are  chemically 
incorporated  into  the  network  in  the  form  of  peroxyl  linkage  Si-O-O-B,  extra  bridging  O  between  two 
B,  and  others.  Since  diffusion  could  be  a  rare  event,  simply  performing  molecular  dynamics  simulations 
and  tracking  the  mean  squared  displacements  may  not  be  sufficient,  and  energy-landscape  exploration 
techniques  such  as  nudged  elastic  band  (NEB)  calculations  may  be  needed.  These  methods  compute  the 
minimum  energy  path  (MEP)  and  saddle-point  configuration  of  thermally  activated  processes,  and  then 
use  transition-state  theory  to  estimate  the  rates. 

First,  we  place  one  Oj  inside  a  cage  and  perform  DFT  molecular  dynamics  simulation  at  2500°C  for  11 
picoseconds.  It  is  very  clear  from  the  MD  trajectory  (in  contrast  to  those  of  simple  liquids  such  as  Ar, 
as  well  as  water)  that  the  borosilicate  liquid  still  maintains  a  very  ’’rigid”  framework  at  2500°C.  Indeed, 
pure  silica  is  a  strong  glass- forming  liquid,  and  its  viscosity  does  not  show  a  precipitous  drop  above  the 
glass  transition  temperature.  For  instance,  even  at  2500°C,  pure  silica  still  has  a  shear  viscosity  ~  104 
Poise,  which  is  a  million  times  thicker  than  that  of  room-temperature  water.  This  means  the  liquid  still 
has  a  well-defined  network  structure  at  any  given  moment,  and  its  topological  change  does  not  occur  at 
the  same  timescale  as,  for  instance,  its  own  Si-O-Si  bond  stretching.  Also,  from  the  11  picoseconds  ab 
initio  MD  trajectory,  we  find  the  02  is  essentially  trapped  inside  one  cage.  It  just  bounces  back  and 
forth  many  times  inside  a  jiggling  cage,  with  no  possibility  of  escape  within  the  MD  simulation  timescale. 
These  facts  suggest  that  one  is  still  justified  to  use  transition-state  theory  and  numerical  schemes  like 
the  nudged  elastic  band  (NEB)  method  to  characterize  diffusion  of  02,  despite  it  is  embedded  in  a 
liquid.  From  our  MD  simulations,  adding  50%  B2O3  to  silica  does  not  seem  to  change  this  consideration 
qualitatively. 

A  complication  for  the  NEB  calculation  is  that  unlike  in  crystals,  diffusion  inside  an  amorphous  frame¬ 
work  has  a  distribution  of  local  minima  and  activation  energies.  Not  all  cages  have  the  same  volume, 
nor  the  same  energy  for  opening  up  the  constrictions  when  OJ  squeezes  from  one  cage  into  the  other. 
In  Fig.  5  of  2,  we  have  shown  a  typical  DFT-NEB  calculated  minimum  energy  path  of  02  diffusion  (ve¬ 
hicular  diffusion  mode).  The  forward  hop  barrier  is  1.8  eV,  whereas  the  backward  hop  barrier  is  1.4  eV. 
These  are  somewhat  higher  than  the  1.12  eV  effective  migration  barrier  that  Bongiorno  and  Pasquarello 
predicted  for  02  vehicular  diffusion  in  pure  silica,  perhaps  because  of  the  B2O3  modifications  to  the 
network.  More  calculations  are  needed  in  order  to  have  better  statistics. 

For  vehicular  diffusion  inside  a  liquid,  there  is  a  well-known  Stokes-Einstein  relation  where  one  correlates 
diffusivity  to  the  inverse  of  viscosity.  Even  though  the  Stokes-Einstein  relation  is  quite  successful  in 
simple  liquids,  it  can  fail  in  network-forming  liquids.  Norton  measured  the  permeation  of  gaseous  oxygen 
through  vitreous  silica  and  found  an  activation  energy  of  27  keal/mol  (1.17  eV)  for  D(0£),  in  substantial 
agreement  with  later  measurements.  However,  the  activation  energy  governing  the  viscosity  of  vitreous 
silica,  is  in  the  range  of  5.3  to  7.5  eV.  So  clearly  the  Stokes-Einstein  relation  does  not  work  here.  From 
our  DFT  modeling,  the  physics  governing  the  activation  energy  of  02  vehicular  diffusion  is  seen  to  be  an 
elastic  deformation  of  the  framework  (elastic  opening  of  the  constrictions)  without  changing  its  network 
topology.  But  the  physics  behind  the  activation  energy  of  must  involve  network  topology  changes,  which 
necessarily  involve  Si-0  bond  breaking. 

2A.  Bongiorno,  C.  J.  Forst,  R.  K.  Kalia,  J.  Li,  J.  Marschall,  A.  Nakano,  M.  M.  Opeka,  I.  G.  Talmy,  P.  Vashishta,  and 
S.  Yip,  MRS  Bulletin,  31  (2006)  410-418. 
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In  addition  to  Oj  vehicular  diffusion,  oxygen  transport  may  also  occur  by  Grotthuss-type  oxygen-hopping 
mechanisms,  mediated  by  network  defects  such  as  peroxyl  linkage  or  oxygen  deficient  centers.  These 
mechanisms  would  involve  bond  breaking,  and  typically  higher  effective  activation  energies  -  mostly  due 
to  the  formation  energies  of  such  defects.  One  such  pathway  involves  an  0£  breaking  up  into  two  O*:  OJ 
— >  0*  +  O*  The  two  0*s  then  move  independently  of  each  other  for  a  while,  and  eventually  recombine 
on  the  other  side  of  the  cage  in  this  particular  NEB  calculation  setup. 

We  find  from  multiple  NEB  calculations  in  the  borosilicate  framework  that  these  coordination  defects, 
once  formed,  can  interconvert  easily,  suggesting  low  migration  barriers,  similar  to  interstitial  defects 
in  metals.  Furthermore,  there  are  two  free  translational  centers:  the  two  CPs,  once  formed,  can  move 
independently  of  each  other  inside  the  liquid.  The  OJ  vehicular  diffusion  has  only  one  free  translational 
center,  in  order  to  maintain  its  molecular  form.  Thus  diffusion  by  OJ  — ►  CP  4-  CP  has  entropic  advantage 
of  approximately  k&  In  c,  where  c  is  the  prevalent  oxygen  carrier  concentration  (per  formula  unit  of 
SiC>2  +  B2O3),  whereas  the  0£  vehicular  diffusion  has  energy  advantage.  The  classic  enthalpy-entropy 
tradeoff  in  free  energy  then  suggests  that  there  exists  a  temperature  7c,  below  which  OJ  is  dominant  in 
concentration,  and  above  which  CPs  are  dominant,  if  the  borosilicate  is  in  an  oxygen-rich  environment 
(the  equivalent  Po2  is  high).  In  lower  equivalent  Pq2  environment,  oxygen  deficient  centers  are  also 
possible  carriers. 

In  reference  to  O2  in  the  gas  phase,  we  find  the  average  potential  energy  of  OJ  solvated  in  borosilicate  is 
0.73  eV  per  molecule,  which  is  essentially  the  elastic  energy  in  the  framework  needed  to  accommodate 
the  molecule.  The  average  potential  energy  of  O*,  on  the  other  hand,  is  about  1.78  eV  per  O.  These 
energies  are  used  to  compute  the  concentration  of  OJ  and  CP  in  borosilicate  in  equilibration  with  Pq2  = 
0.1  atm.  Note  that  0(0^)  scales  linearly  with  Pq2,  whereas  c(CP)  scales  linearly  with  Pq2  ■  Thus,  deeper 
and  deeper  into  the  scale,  as  the  oxygen  chemical  potential  gets  lower  and  lower,  the  Grotthuss-type 
oxygen  transport  should  become  relatively  more  and  more  important. 

Norton  measured  the  solubility  of  in  pure  silica  at  1078°C  and  Po2  =  1  atm  to  be  1.9  x  10-3  cm3 
STP  O2  gas  per  cm3  silica.  This  amounts  to  5.1  x  1016  per  cm3  silica  or  a  dimensionless  concentration 
of  about  c  =  10-5,  at  Pq2  =  1  atm.  This  is  about  two  orders  of  magnitude  higher  than  our  DFT 
predicted  0£  solubility  in  borosilicate.  This  could  be  due  to  the  structural  difference  between  pure  silica 
framework  and  borosilicate  framework.  As  Bongiorno  and  Pasquarello  noted,  the  potential  energy  of  Oj 
depends  sensitively  on  the  cage  interstice  volume,  and  there  is  certainly  a  structural  difference  between 
pure  silica  and  B2O3  modified  networks.  This  could  also  be  partly  due  to  the  intrinsic  errors  of  PW91 
density  functional,  which  are  known  to  give  large  errors  treating  isolated  molecules  (the  reference  state), 
and  non-bonding  interactions  (Oj  interactions  with  the  framework). 

The  above  uncertainties  aside,  it  is  still  a  rather  conservative  estimate  that  above  Tc  ~  1500°C  the 
dominant  oxygen  carriers  in  the  borosilicate  "liquid  roof’  are  network  defects  instead  of  molecular  OJ, 
because  there  are  other  factors  that  disfavors  the  Deal-Grove  mechanism:  (a)  the  DFT  calculations 
indicate  that  the  CPs  not  only  have  entropy  advantage,  but  also  mobility  advantages  over  OJ,  (b) 
intense  aerothermal  heating  environment  may  introduce  a  significant  level  of  dissociated  oxygen  on  the 
outer  liquid  surface,  which  would  favor  O*  diffusion  from  a  non-equilibrium  kinetics  perspective,  (c)  as 
Po2  drops  from  ~  104  Pa  on  the  outer  liquid  surface  to  less  than  10-5  Pa  equivalent  at  the  internal  gas 
finger-liquid  interface  according  to  the  volatility  diagram,  the  balance  will  shift  more  and  more  away 
from  OJ  vehicular  diffusion  to  network  defect  Grotthuss  diffusion.  Oxygen  deficient  centers  may  play 
a  significant  role  in  oxygen  transport  at  low  equivalent  Po2s  ^  seems  plausible  that  oxygen  in  the 
borosilicate  liquid  could  react  with  the  underlying  substrate  or  SiO  (g),  injecting  oxygen  vacancies  (e.g. 
regions  of  high  B  and  Si  stoichiometry)  into  the  liquid,  which  then  diffuse  up  the  scale  to  recombine 
with  Oj  or  CP  somewhere  inside  the  liquid,  (d)  In  all  the  DFT  calculations,  we  have  only  considered 
neutral  network  defects.  Charging  the  defects  may  significantly  greatly  reduce  their  formation  energies, 
although  we  will  then  need  to  solve  the  complementary  problem  of  what  other  defects  compensate  the 
charge  and  carry  out  arnbipolar  diffusion  under  open-circuit  condition. 
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The  mobility  advantage  of  the  peroxyl  linkage  over  Oj  suggests  that  arc-jet  testing,  which  introduces 
a  non-equilibrium  distribution  of  dissociated  oxygen  atoms  on  the  surface,  will  likely  lead  to  faster 
oxidation  than  ordinary  furnace  testing  at  the  same  temperature.  Also,  because  the  network  defects 
are  chemically  incorporated  into  the  network  and  thus  interact  more  strongly  with  solutes  than  OJ, 
small  changes  in  the  glass  chemistry  could  lead  to  large  changes  in  the  oxygen  diffusivity  by  defect 
trapping/gettering. 

Summary 

We  present  a  congruent  explanation  of  the  oxidation  protection  of  ZrB2+SiC  based  on  a  “solid  pillars, 
liquid  roof”  scale  architecture,  where  the  borosilicate  liquid  defines  the  effective  diffusion  barrier,  and  the 
solid  zirconia  collects  and  retains  the  liquid  and  provides  mechanical  support.  Internal  gas  fingers  will 
form  as  the  liquid  phase  retreats  to  remove  the  high  SiO  volatility  above  a  boiling  transition  temperature. 

At  such  high  temperatures,  to  satisfy  both  thermochemical  and  mechanical  stabilities,  the  “solid  pillars, 
liquid  roof”  architecture  seems  to  be  a  viable  solution,  not  available  to  monolithic  SiC. 

Compared  to  the  borosilicate  liquid  phase,  whether  fully  dense  zirconia  is  blocking  or  unblocking  to 
oxygen  in  open-circuit  condition  depends  on  its  electronic  conductivity  (transference  number),  which  in 
turn  depends  on  how  the  charge  defects  are  compensated  inside  the  crystal,  related  to  the  amount  of 
impurities.  If  the  zirconia  phase  has  a  highly  porous  microstructure,  however,  then  the  above  discussion 
is  likely  irrelevant  and  the  borosilicate  liquid  phase  will  control  the  effective  diffusion  barrier,  since  it 
will  occupy  both  serial  and  parallel  oxygen  transport  pathways. 

At  low  temperatures,  it  is  commonly  accepted  that  molecular  oxygen  OJ  dominates  oxygen  transport. 
However,  from  first-principles  calculations  with  detailed  borosilicate  atomic  structures,  it  seems  unlikely 
that  this  will  remain  the  case  at  temperatures  of  practical  interest  for  the  ZrB2+SiC  thermal  protec¬ 
tion  system  (above  1500°C).  This  means  that  the  oxidation  rate  will  likely  have  a  complex,  sub-linear 
dependence  with  respect  to  the  external  oxygen  partial  pressure.  Also,  if  the  oxygen  carriers  are  chem¬ 
ically  incorporated  and  interact  strongly  with  the  framework,  there  is  hope  that  by  tuning  the  glass 
composition,  the  carriers  could  be  trapped,  thereby  slowing  down  oxygen  diffusion. 
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